arXiv:1508.07453vl [astro-ph.HE] 29 Aug 2015 


PROCEEDINGS 

OF SCIENCE 



A detailed study of Gamma-ray emission from 
Cassiopeia A using VERITAS 


Sajan Kumar^ for VERITAS collaboration* *^ 

^University of Delaware, 19716, USA 
E-mail: sajan@udel.edu 

Supernova remnants (SNRs) have long been considered the leading candidate sites for the ac¬ 
celeration of cosmic rays within the Galaxy through the process of diffusive shock acceleration. 
The connection between SNRs and cosmic rays is supported by the detection of high energy (HE; 
100 MeV to 100 GeV) and very high energy (VHE; 100 GeV to 100 TeV) gamma rays from 
young and middle-aged SNRs. However, the interpretation of the gamma-ray observations is not 
unique. This is due to the fact that gamma rays can be produced by electrons through non-thermal 
Bremsstrahlung and inverse Compton scattering, and by protons through proton-proton collisions 
and subsequent neutral pion decay. To disentangle and quantify the contributions of electrons and 
protons to the gamma-ray flux, it is necessary to measure precisely the spectra and morphology of 
SNRs over a broad range of gamma-ray energies. Cassiopeia A (Cas A) is one such young SNR 
( 350 years) which is bright in radio and X-rays. It has been detected as a bright point source in 
HE gamma rays by Eermi-LAT and in VHE gamma rays by HEGRA, MAGIC and VERITAS. 
Cas A has been observed with VERITAS for more than 60 hours - three times the published ex¬ 
posure. The observations were taken between 2007 and 2013 over a wide range of zenith angles. 
In particular, half of the total data were taken under large zenith angles to boost the effective area 
above few TeV. Here, we will present the detailed spectral and morphological results from the 
complete dataset. 
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1. Introduction 

For decades, young supernova remnants(SNRs) have been thought to be responsible for the 
production of the bulk of Galactic cosmic rays through the mechanism of diffusive shock accel¬ 
eration. The support for this conviction is backed by a number of observational facts such as 
the detection of non-thermal X-rays and high energy/very high energy gamma rays from these 
SNRs[l]. In the class of SNRs, Cassiopeia A (Cas A) is one of the youngest (age ~ 350 yrs [2]) 
object, located at a distance of (~ 3.4 kpc [3]). Observations of the scattered light echo from the 
supernova explosion indicate that Cas A was a type Ilb supernova from a star of initial mass 15-25 
Mq and had lost most of its hydrogen envelop prior to exploding [4] . Due to its proximity to Earth 
and its well established progenitor and brightness, it has been observed extensively across the elec¬ 
tromagnetic spectrum, from radio to gamma rays. Viewed in radio wavelengths [5,6], it is one of 
the brightest sources, with the main fraction of the emission coming from a bright radio ring (radius 
~ 100") thought to be associated with the reverse shock. This radio emission is generated through 
the synchrotron process by electrons that are trapped and accelerated in the magnetic field. Optical 
emission is dominated mainly by dense ejecta heated by reverse-shock and several fast moving 
knots [7]. Broadening of CO line emission towards Cas A has been observed by the Spitzer space 
telescope, which suggests an interaction between the shock front of the SNR and molecular clouds 
[ 8 , 9]. In X-rays (~ 0.1 keV to 80 keV), thermal and non-thermal emission is seen with XMM- 
Newton, Chandra and NuSTAR. Thermal X-ray emission arises from the shell of reverse-shocked 
ejecta, rich in emission lines from highly ionized atoms and overlapping with the bright radio shell 
[10, 11]. Besides this, non-thermal emission is also seen at both forward and reverse shock [12, 13] 
, although the emission at the reverse shock can also be interpreted as the emission from the for¬ 
ward shock projected towards the inner ring [14]. Recently, observations from NuSTAR resolve 
the remnant above 15 KeV and find fhaf emission is dominafed nof only by forward and reverse 
shocks buf also by knofs locafed in fhe inferior of remnanf [15]. This non-fhermal X-ray emission 
can be inferprefed as synchrofron radiafion from elecfrons accelerafed fo a few fens of TeV wifhin 
shocks. 

Since elecfrons are being accelerafed fo relafivisfic energies af bofh forward and reverse shocks, 
if has been argued fhaf gamma-rays (from MeV fo TeV range) should be expecfed from fhis source 
fhrough inverse Compfon scaffering. However, if is difficulf fo differenfiafe fhe confribufion of fhe 
forward and reverse shocks fo gamma-ray emission because of fhe limifed angular resolution of 
space-based and ground-based gamma-ray observafories. Alfernafively, if profons are accelerafed 
af fhe shock fronfs, fhey can inferacf wifh ambienf fargef maferial (gas), also producing gamma-ray 
emission fhrough neufral pion decay. Recenfly, direcf evidence for fhe accelerafion of profons af 
SNR shocks has been confirmed from fwo SNRs: IC443 and W44 [16]. 

In fhe HE gamma-ray energy range (MeV fo GeV), fhe Fermi Gamma ray space Telescope 
defected a signal from Cas A for fhe firsf fime in 2010 [17], buf did nof provide very sfrong con- 
sfrainfs on fhe emission models. A more recenf paper by Yuan el al. [18] shows fhaf a hadronic 
model is preferred over leplonic model in fhe GeV energy range and fhis poinl is supported by Saha 
el al. [19] also. In TeV regime, HEGRA firsf discovered VHE gamma-ray emission from Cas A 
[20], wifh laler confirmalions of fhis deleclion supplied by MAGIC and VERITAS observatories 
[21, 22]. The speclral index (7 = 2.4 — 2.6) and fhe fluxes( 3% of fhe Crab Nebula flux) published 
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by these three groups match well with each other within errors. Due to the limited angular reso¬ 
lution of these instruments, Cas A is detected as a point source of gamma rays in the TeV energy 
range. 

In this work, we will present the results from observations of Cas A with VERITAS taken 
between 2007 and 2013, which amount to more than 60 hours. This exposure is almost three 
times that already published by VERITAS, and significantly reduces the statistical error on the 
flux and spectral index, providing further constraints which help to distinguish between hadronic 
and leptonic models of gamma-ray emission. Here we present a brief summary of the VERITAS 
instrument, followed by the procedures for data selection and analysis, and discuss the results and 
conclusions. 


2. VERITAS Telescope 

VERITAS (the Very Energetic Radiation Imaging Telescope Array System) is a ground-based 
gamma-ray observatory which consists of an array of four telescopes, located in southern arizona at 
(30°40'30.21"N, 110°57'7.77"W). Each telescope has a 12-m diameter optical reflector, providing 
a total reflecting area of ~ llOm^. A camera consisting of 499 photomultiplier tubes (PMTs) is 
placed in the focal plane of each telescope, giving a field of view of 3.5°. VERITAS employs the 
technique of stereoscopic reconstruction of gamma-ray initiated air showers. Erom 2007 to 2012, 
the array underwent two major upgrades. The first occurred during the summer of 2009, when the 
original prototype telescope was moved to a new location to make the whole array layout more 
symmetric and improve the baseline distances. In the second upgrade, in 20II-20I2, EPGA-based 
E2 trigger system and higher efficiency PMTs (super-bialkali photocathod) were installed [23, 24] . 
Along with the inclusion of new data analysis techniques, these upgrades have resulted in improved 
sensitivity of the array and a lower energy threshold. Currently, a source with a flux level of 1 % of 
the Crab Nebula can be detected in 25h. The angular resolution for gamma-rays at 1 TeV energy is 
~ 0.08° and the sensitive energy detection range spans from 85 GeV to > 30 TeV. 


3. Data selection and Analysis 

The observations of Cas A are summarised in Table 1. All observations were performed using 
four telescopes under very dark and clear sky conditions. After quality selection cuts the total ex¬ 
posure is ~65 hours. The first step in the analysis procedure is image calibration and cleaning [25]. 
This step is necessary to select only those pixels which contain the Cherenkov light and remove all 
other pixels containing night sky background. Calculation of Hillas parameters (length, width and 
size of image) is the second step [26]. Prom these parameters, a gamma ray initiated shower can be 
distinguished from a background shower initiated by cosmic rays. After the image parametrization, 
the aiTival direction of each event is reconstructed using the stereoscopic technique, wherein the 
major axes of the shower images intersect at a point in the camera plane, providing a geometrical 
method to measure the arrival direction of an event. 
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Table 1: Details of observations of Cas A. 


Data Set 

Date 


dz range 
(deg) 

Average dz 
(deg) 

Wobble 

(deg) 

Eive Time 

(Hours) 

Mean trigger rate 
(Hz) 

I 

09/07 - 11/07 

4 

27-40 

34 

0.5 

18 

250 

II 

12/11 - 12/11 

4 

33-43 

38 

0.5 

2 

350 

III 

09/12 - 12/13 

4 

24-39 

30 

0.5 

19 

400 

IV 

09/12 - 12/13 

4 

40-64 

56 

0.5 

25 

300 


4. Results 

The VERITAS skymap of excess VHE gamma-ray events from the direction of Cas A is 
shown in Eigure 1 (left). This skymap is generated using the reflected background model [27] 
and smoothed with a circular window of radius 0.09°. Note that the skymap shown here results 
from the analysis of 2012-2013 data (taken at small zenith angles) and corresponds to a statisti¬ 
cal significance of 11 a. The position of the gamma-ray source is determined by fitting a two- 
dimensional Gaussian function to a non-smoothed skymap. The blue cross on the skymap with 
RA=23''23'"20.4"±0°.006 ,m ±0°.01 V and Dec= 58.817 ±0°.006,m ±0°.014,y, shows the cen¬ 
troid of the gamma-ray source. Erom this analysis, it is clear that position of the gamma-ray source 
is mainly limited by the systematic error in the pointing of the telescopes, even for this reduced 
dataset. In Eigure 1 (right), we shows the comparison of the centroid locations from the Fermi 
(Yuan et al 2013 [18]), VERITAS (Acciari et al 2010 [22]), MAGIC (Albert et al 2007 [21]) and 
VERITAS (new work). The new VERITAS centroid for the TeV emission is consistent with the 
center of the remnant, with the other TeV measurements, and with the Eermi-EAT position. The 
radius of GeV and TeV regions corresponds to 1 a statistical and systematic errors added in quadra¬ 
ture. 

The differential energy spectrum for the whole data set measured by VERITAS is also shown 
in Eigure 2, along with the published data points from VERITAS. The spectral points are fitted with 
a power-law in the energy range from 300 GeV to 7 TeV: 

^ = (1.45±0.11) X 

Cl 12/ 

A power-law fit to the data points gives a of 2.22 for 5 degrees of freedom, resulting in a 
good fit probability of 81%. 


5. Summary and Discussions 

Cas A has been observed by VERITAS for ~ 65 hours over a period of 6 years. One third of 
the total data has already been published [22]. By adding the new data to the already published 
data, we refine the spectral measurements (from our published spectrum) both at lower and higher 
energy. In addition to this, the improved signal-to-noise ratio helps to reduce the statistical errors 
on the spectral index by more than 50%. Eor the centroid of the gamma-ray emission, we are 
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Figure 1: (left) The skymap of excess counts from the region of Cas A, smoothed with a circular window 
of radius 0.09°. This map was produced using 18 hours of VERITAS observations from 2012 (with the 
upgraded camera and at small zenith angles). The white circle indicates the size of VERITAS point spread 
function. The cross indicates the measured position of the TeV gamma-ray source, (right) Comparison of 
centroid positions from Fermi (yellow, [18]), VERITAS (green, [22]) and MAGIC (red, [21]) with the new 
VERITAS (white) centroid position. 

now limited by the systematics in the pointing of our telescopes (50 arcseconds). Figure 3 shows 
the complete gamma-ray spectral energy distribution using the new flux points from VERITAS and 
those from Fermi [18]. A discussion of these results will be presented at the conference. 
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Figure 2: Differential energy spectrum from Cas A (this work) along with published VERITAS points (blue, 
Acciari et al 2010 [22]). 
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Figure 3: Combined differential flux points using Fermi (Yuan et al 2013 [18]) and VERITAS data (this 
work). 
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